Using analytical formulation of Bogner-Fox-Schmit (BFS) plate element
Span lengths
โa = hpโ(โ[3.6; 4.2; 4.2; 3.6]โ)โ = [3.6โ4.2โ4.2โ3.6] m
โb = hpโ(โ[3; 3.6; 3]โ)โ = [3โ3.6โ3] m
Number of axes - nsa = lenโ(โโaโ)โ + 1 = 5 , nsb = lenโ(โโbโ)โ + 1 = 4
Axis coordinates - โxs = [0โ3.6โ7.8โ12โ15.6] m, โys = [0โ3โ6.6โ9.6] m
Slab dimensions - la = โxs.5 = 15.6 m, lb = โys.4 = 9.6 m
Thickness - t = 0.2 m
Load - q = 10 kN/mยฒ
Modulus of elasticity - E = 35000 MPa
Poisson`s ratio - ฮฝ = 0.2
We will use BFS rectangular finite element with nDOFs = 16
Element dimensions - a1 = 0.6 m, b1 = 0.6 m
Number of elements and joints along a and b -
โna = ceiling(โaa1) = ceiling(โa0.6) = [6โ7โ7โ6] , nea = sumโ(โโnaโ)โ = 26 , nja = nea + 1 = 26 + 1 = 27
โnb = ceiling(โbb1) = ceiling(โb0.6) = [5โ6โ5] , neb = sumโ(โโnbโ)โ = 16 , njb = neb + 1 = 16 + 1 = 17
Total number of elements - ne = neaโยทโneb = 26โยทโ16 = 416
Total number of joints - nj = njaโยทโnjb = 27โยทโ17 = 459
Supported joints count - ns = nsaโยทโnsb = 5โยทโ4 = 20
Joint coordinates
โxj = [0โ0โ0โ0โ0โ0โ0โ0โ0โ0โ0โ0โ0โ0โ0โ0โ0โ0.6โ0.6โ0.6โ...โ15.6] m
โyj = [0โ0.6โ1.2โ1.8โ2.4โ3โ3.6โ4.2โ4.8โ5.4โ6โ6.6โ7.2โ7.8โ8.4โ9โ9.6โ0โ0.6โ1.2โ...โ9.6] m
Numbers of joints at elements' corners
transpโ(โejโ)โ = 1234567891011121314151618192021โฏ441 1819202122232425262728293031323335363738โฏ458 1920212223242526272829303132333436373839โฏ459 23456789101112131415161719202122โฏ442
Supported joints
โsj = [1โ6โ12โ17โ103โ108โ114โ119โ222โ227โ233โ238โ341โ346โ352โ357โ443โ448โ454โ459]
Joints for element e - jeโ(โeโ)โ = rowโ(โej; eโ)โ
Shape functions
Along dimension a
| Base functions
ฮฆ1aโ(โฮพโ)โ = 1 โ ฮพ2โยทโโ(โ3 โ 2โยทโฮพโ)โ ฮฆ2aโ(โฮพโ)โ = ฮพโยทโa1โยทโโ(โ1 โ ฮพโยทโโ(โ2 โ ฮพโ)โโ)โ ฮฆ3aโ(โฮพโ)โ = ฮพ2โยทโโ(โ3 โ 2โยทโฮพโ)โ ฮฆ4aโ(โฮพโ)โ = ฮพ2โยทโa1โยทโโ(โ-1 + ฮพโ)โ | First derivatives
ฮฆโฒ1aโ(โฮพโ)โ = -6โยทโฮพa1โยทโโ(โ1 โ ฮพโ)โ ฮฆโฒ2aโ(โฮพโ)โ = 1 โ ฮพโยทโโ(โ4 โ 3โยทโฮพโ)โ ฮฆโฒ3aโ(โฮพโ)โ = 6โยทโฮพa1โยทโโ(โ1 โ ฮพโ)โ ฮฆโฒ4aโ(โฮพโ)โ = -ฮพโยทโโ(โ2 โ 3โยทโฮพโ)โ | Second derivatives
ฮฆโณ1aโ(โฮพโ)โ = -โ6a12โยทโโ(โ1 โ 2โยทโฮพโ)โ ฮฆโณ2aโ(โฮพโ)โ = -โ2a1โยทโโ(โ2 โ 3โยทโฮพโ)โ ฮฆโณ3aโ(โฮพโ)โ = 6a12โยทโโ(โ1 โ 2โยทโฮพโ)โ ฮฆโณ4aโ(โฮพโ)โ = -โ2a1โยทโโ(โ1 โ 3โยทโฮพโ)โ |
Along dimension b
| Base functions
ฮฆ1bโ(โฮทโ)โ = 1 โ ฮท2โยทโโ(โ3 โ 2โยทโฮทโ)โ ฮฆ2bโ(โฮทโ)โ = ฮทโยทโb1โยทโโ(โ1 โ ฮทโยทโโ(โ2 โ ฮทโ)โโ)โ ฮฆ3bโ(โฮทโ)โ = ฮท2โยทโโ(โ3 โ 2โยทโฮทโ)โ ฮฆ4bโ(โฮทโ)โ = ฮท2โยทโb1โยทโโ(โ-1 + ฮทโ)โ | First derivatives
ฮฆโฒ1bโ(โฮทโ)โ = -6โยทโฮทb1โยทโโ(โ1 โ ฮทโ)โ ฮฆโฒ2bโ(โฮทโ)โ = 1 โ ฮทโยทโโ(โ4 โ 3โยทโฮทโ)โ ฮฆโฒ3bโ(โฮทโ)โ = 6โยทโฮทb1โยทโโ(โ1 โ ฮทโ)โ ฮฆโฒ4bโ(โฮทโ)โ = -ฮทโยทโโ(โ2 โ 3โยทโฮทโ)โ | Second derivatives
ฮฆโณ1bโ(โฮทโ)โ = -โ6b12โยทโโ(โ1 โ 2โยทโฮทโ)โ ฮฆโณ2bโ(โฮทโ)โ = -โ2b1โยทโโ(โ2 โ 3โยทโฮทโ)โ ฮฆโณ3bโ(โฮทโ)โ = 6b12โยทโโ(โ1 โ 2โยทโฮทโ)โ ฮฆโณ4bโ(โฮทโ)โ = -โ2b1โยทโโ(โ1 โ 3โยทโฮทโ)โ |
| For vertical displacements w
N1,wโ(โฮพ; ฮทโ)โ = ฮฆ1aโโ(โฮพโ)โโยทโฮฆ1bโโ(โฮทโ)โ N2,wโ(โฮพ; ฮทโ)โ = ฮฆ3aโโ(โฮพโ)โโยทโฮฆ1bโโ(โฮทโ)โ N3,wโ(โฮพ; ฮทโ)โ = ฮฆ3aโโ(โฮพโ)โโยทโฮฆ3bโโ(โฮทโ)โ N4,wโ(โฮพ; ฮทโ)โ = ฮฆ1aโโ(โฮพโ)โโยทโฮฆ3bโโ(โฮทโ)โ | For rotations ฮธโ
N1,ฮธโโ(โฮพ; ฮทโ)โ = ฮฆ2aโโ(โฮพโ)โโยทโฮฆ1bโโ(โฮทโ)โ N2,ฮธโโ(โฮพ; ฮทโ)โ = ฮฆ4aโโ(โฮพโ)โโยทโฮฆ1bโโ(โฮทโ)โ N3,ฮธโโ(โฮพ; ฮทโ)โ = ฮฆ4aโโ(โฮพโ)โโยทโฮฆ3bโโ(โฮทโ)โ N4,ฮธโโ(โฮพ; ฮทโ)โ = ฮฆ2aโโ(โฮพโ)โโยทโฮฆ3bโโ(โฮทโ)โ | For rotations ฮธแตง
N1,ฮธแตงโ(โฮพ; ฮทโ)โ = ฮฆ1aโโ(โฮพโ)โโยทโฮฆ2bโโ(โฮทโ)โ N2,ฮธแตงโ(โฮพ; ฮทโ)โ = ฮฆ3aโโ(โฮพโ)โโยทโฮฆ2bโโ(โฮทโ)โ N3,ฮธแตงโ(โฮพ; ฮทโ)โ = ฮฆ3aโโ(โฮพโ)โโยทโฮฆ4bโโ(โฮทโ)โ N4,ฮธแตงโ(โฮพ; ฮทโ)โ = ฮฆ1aโโ(โฮพโ)โโยทโฮฆ4bโโ(โฮทโ)โ |
|
For twist ฯ N1,ฯโ(โฮพ; ฮทโ)โ = ฮฆ2aโโ(โฮพโ)โโยทโฮฆ2bโโ(โฮทโ)โ N2,ฯโ(โฮพ; ฮทโ)โ = ฮฆ4aโโ(โฮพโ)โโยทโฮฆ2bโโ(โฮทโ)โ N3,ฯโ(โฮพ; ฮทโ)โ = ฮฆ4aโโ(โฮพโ)โโยทโฮฆ4bโโ(โฮทโ)โ N4,ฯโ(โฮพ; ฮทโ)โ = ฮฆ2aโโ(โฮพโ)โโยทโฮฆ4bโโ(โฮทโ)โ |
| |
|
N1,w shape function plot |
N1,ฮธโ shape function plot |
Constitutive matrix (stress - strain relationship)
D1,1 = Eโยทโt312โยทโโ(โ1 โ ฮฝ2โ)โ = 35000โยทโ0.2312โยทโโ(โ1 โ 0.22โ)โ = 24.305556 kNm
D = D1,1โยทโhp([1; ฮฝ; 0 | ฮฝ; 1; 0 | 0; 0; 1 โ ฮฝ2]) = 24.305556โยทโhp([1; 0.2; 0 | 0.2; 1; 0 | 0; 0; 1 โ 0.22]) = 24.3055564.8611110 4.86111124.3055560 009.722222 kNm
ฮฑ = b1a1 = 0.60.6 = 1
Ke.1, 1 = 15635โยทโ(ฮฑa12 + 1ฮฑโยทโb12) + 7225โยทโA1 = 15635โยทโ(10.62 + 11โยทโ0.62) + 7225โยทโ0.36 = 32.761905
Ke.5, 5 = Ke.1, 1 = 32.761905 , Ke.9, 9 = Ke.1, 1 = 32.761905 , Ke.13, 13 = Ke.1, 1 = 32.761905
Ke.1, 2 = 235โยทโ(39โยทโฮฑa1 + 11ฮฑยฒโยทโb1) + 30โยทโฮฝ + 625โยทโb1 = 235โยทโ(39โยทโ10.6 + 111โยทโ0.6) + 30โยทโ0.2 + 625โยทโ0.6 = 5.561905
Ke.5, 6 = -Ke.1, 2 = -5.561905 , Ke.9, 10 = -Ke.1, 2 = -5.561905 , Ke.13, 14 = Ke.1, 2 = 5.561905
Ke.1, 3 = 235โยทโ(11โยทโฮฑยฒa1 + 39ฮฑโยทโb1) + 30โยทโฮฝ + 625โยทโa1 = 235โยทโ(11โยทโ10.6 + 391โยทโ0.6) + 30โยทโ0.2 + 625โยทโ0.6 = 5.561905
Ke.5, 7 = Ke.1, 3 = 5.561905 , Ke.9, 11 = -Ke.1, 3 = -5.561905 , Ke.13, 15 = -Ke.1, 3 = -5.561905
Ke.1, 4 = 1135โยทโ(ฮฑยฒ + 1ฮฑยฒ) + 10โยทโฮฝ + 150 = 1135โยทโ(1 + 11) + 10โยทโ0.2 + 150 = 0.6885714
Ke.5, 8 = -Ke.1, 4 = -0.6885714 , Ke.9, 12 = Ke.1, 4 = 0.6885714 , Ke.13, 16 = -Ke.1, 4 = -0.6885714
Ke.1, 5 = 235โยทโ(27ฮฑโยทโb12 โ 78โยทโฮฑa12) โ 7225โยทโA1 = 235โยทโ(271โยทโ0.62 โ 78โยทโ10.62) โ 7225โยทโ0.36 = -16.095238
Ke.9, 13 = Ke.1, 5 = -16.095238
Ke.1, 6 = 1335โยทโ(6โยทโฮฑa1 โ 1ฮฑยฒโยทโb1) + 625โยทโb1 = 1335โยทโ(6โยทโ10.6 โ 11โยทโ0.6) + 625โยทโ0.6 = 3.495238
Ke.2, 5 = -Ke.1, 6 = -3.495238 , Ke.9, 14 = -Ke.1, 6 = -3.495238 , Ke.10, 13 = Ke.1, 6 = 3.495238
Ke.1, 7 = 135โยทโ(27ฮฑโยทโb1 โ 22โยทโฮฑยฒa1) โ 30โยทโฮฝ + 625โยทโa1 = 135โยทโ(271โยทโ0.6 โ 22โยทโ10.6) โ 30โยทโ0.2 + 625โยทโ0.6 = -0.5619048
Ke.3, 5 = Ke.1, 7 = -0.5619048 , Ke.9, 15 = -Ke.1, 7 = 0.5619048 , Ke.11, 13 = -Ke.1, 7 = 0.5619048
Ke.2, 4 = 235โยทโ(11โยทโฮฑโยทโb13 + a1ฮฑยฒ) + 2โยทโa1โยทโโ(โ5โยทโฮฝ + 1โ)โ75 = 235โยทโ(11โยทโ1โยทโ0.63 + 0.61) + 2โยทโ0.6โยทโโ(โ5โยทโ0.2 + 1โ)โ75 = 0.192
Ke.6, 8 = Ke.2, 4 = 0.192 , Ke.10, 12 = -Ke.2, 4 = -0.192 , Ke.14, 16 = -Ke.2, 4 = -0.192
Ke.2, 6 = 135โยทโ(26โยทโฮฑ โ 3ฮฑ3) โ 225โยทโฮฑ = 135โยทโ(26โยทโ1 โ 313) โ 225โยทโ1 = 0.5771429
Ke.10, 14 = Ke.2, 6 = 0.5771429
Ke.2, 8 = 135โยทโ(11โยทโฮฑโยทโb13 โ 3โยทโa12โยทโฮฑยฒ) โ a1โยทโโ(โ5โยทโฮฝ + 1โ)โ150 = 135โยทโ(11โยทโ1โยทโ0.63 โ 3โยทโ0.62โยทโ1) โ 0.6โยทโโ(โ5โยทโ0.2 + 1โ)โ150 = 0.02914286
Ke.4, 6 = Ke.2, 8 = 0.02914286 , Ke.10, 16 = -Ke.2, 8 = -0.02914286 , Ke.12, 14 = -Ke.2, 8 = -0.02914286
Ke.2, 10 = 335โยทโ(1ฮฑ3 + 3โยทโฮฑ) + 225โยทโฮฑ = 335โยทโ(113 + 3โยทโ1) + 225โยทโ1 = 0.4228571
Ke.6, 14 = Ke.2, 10 = 0.4228571
Ke.2, 12 = -โ170โยทโ(3โยทโa1ฮฑยฒ + 13โยทโฮฑโยทโb13) โ a1150 = -โ170โยทโ(3โยทโ0.61 + 13โยทโ1โยทโ0.63) โ 0.6150 = -0.06685714
Ke.4, 10 = -Ke.2, 12 = 0.06685714 , Ke.6, 16 = Ke.2, 12 = -0.06685714 , Ke.8, 14 = -Ke.2, 12 = 0.06685714
Ke.2, 14 = 235โยทโ(9โยทโฮฑ โ 2ฮฑ3) โ 825โยทโฮฑ = 235โยทโ(9โยทโ1 โ 213) โ 825โยทโ1 = 0.08
Ke.6, 10 = Ke.2, 14 = 0.08
Ke.2, 15 = 135โยทโ(11ฮฑยฒ โ 13โยทโฮฑยฒ2) + 5โยทโฮฝ + 150 = 135โยทโ(111 โ 13โยทโ12) + 5โยทโ0.2 + 150 = 0.1685714
Ke.3, 14 = -Ke.2, 15 = -0.1685714 , Ke.6, 11 = -Ke.2, 15 = -0.1685714 , Ke.7, 10 = Ke.2, 15 = 0.1685714
Ke.1, 14 = 135โยทโ(27โยทโฮฑa1 โ 22ฮฑยฒโยทโb1) โ 30โยทโฮฝ + 625โยทโb1 = 135โยทโ(27โยทโ10.6 โ 221โยทโ0.6) โ 30โยทโ0.2 + 625โยทโ0.6 = -0.5619048
Ke.2, 13 = Ke.1, 14 = -0.5619048 , Ke.5, 10 = -Ke.1, 14 = 0.5619048 , Ke.6, 9 = -Ke.1, 14 = 0.5619048
Ke.1, 15 = 1335โยทโ(6ฮฑโยทโb1 โ ฮฑยฒa1) + 625โยทโa1 = 1335โยทโ(61โยทโ0.6 โ 10.6) + 625โยทโ0.6 = 3.495238
Ke.3, 13 = -Ke.1, 15 = -3.495238 , Ke.5, 11 = Ke.1, 15 = 3.495238 , Ke.7, 9 = -Ke.1, 15 = -3.495238
Ke.1, 16 = 135โยทโ(11ฮฑยฒ โ 13โยทโฮฑยฒ2) + 5โยทโฮฝ + 150 = 135โยทโ(111 โ 13โยทโ12) + 5โยทโ0.2 + 150 = 0.1685714
Ke.4, 13 = -Ke.1, 16 = -0.1685714 , Ke.8, 9 = Ke.1, 16 = 0.1685714 , Ke.5, 12 = -Ke.1, 16 = -0.1685714
Ke.2, 2 = 435โยทโ(13โยทโฮฑ + 1ฮฑ3) + 825โยทโฮฑ = 435โยทโ(13โยทโ1 + 113) + 825โยทโ1 = 1.92
Ke.6, 6 = Ke.2, 2 = 1.92 , Ke.10, 10 = Ke.2, 2 = 1.92 , Ke.14, 14 = Ke.2, 2 = 1.92
Ke.2, 3 = 1135โยทโ(ฮฑยฒ + 1ฮฑยฒ) + 60โยทโฮฝ + 150 = 1135โยทโ(1 + 11) + 60โยทโ0.2 + 150 = 0.8885714
Ke.6, 7 = -Ke.2, 3 = -0.8885714 , Ke.10, 11 = Ke.2, 3 = 0.8885714 , Ke.14, 15 = -Ke.2, 3 = -0.8885714
Ke.3, 16 = -โ135โยทโ(3โยทโฮฑยฒโยทโb12 โ 11โยทโa13โยทโฮฑ) โ b1โยทโโ(โ5โยทโฮฝ + 1โ)โ150 = -โ135โยทโ(3โยทโ1โยทโ0.62 โ 11โยทโ0.63โยทโ1) โ 0.6โยทโโ(โ5โยทโ0.2 + 1โ)โ150 = 0.02914286
Ke.4, 15 = Ke.3, 16 = 0.02914286 , Ke.7, 12 = -Ke.3, 16 = -0.02914286 , Ke.8, 11 = -Ke.3, 16 = -0.02914286
Ke.4, 4 = 4105โยทโ(ฮฑโยทโb12 + a12ฮฑ) + 8โยทโA1225 = 4105โยทโ(1โยทโ0.62 + 0.621) + 8โยทโ0.36225 = 0.04022857
Ke.8, 8 = Ke.4, 4 = 0.04022857 , Ke.12, 12 = Ke.4, 4 = 0.04022857 , Ke.16, 16 = Ke.4, 4 = 0.04022857
Ke.4, 8 = 135โยทโ(2โยทโฮฑโยทโb123 โ a12ฮฑ) โ 2โยทโA1225 = 135โยทโ(2โยทโ1โยทโ0.623 โ 0.621) โ 2โยทโ0.36225 = -0.006628571
Ke.12, 16 = Ke.4, 8 = -0.006628571
Ke.4, 12 = -โ170โยทโ(a12ฮฑ + ฮฑโยทโb12) + A1450 = -โ170โยทโ(0.621 + 1โยทโ0.62) + 0.36450 = -0.009485714
Ke.8, 16 = Ke.4, 12 = -0.009485714
Ke.4, 16 = 135โยทโ(2โยทโa123โยทโฮฑ โ ฮฑโยทโb12) โ 2โยทโA1225 = 135โยทโ(2โยทโ0.623โยทโ1 โ 1โยทโ0.62) โ 2โยทโ0.36225 = -0.006628571
Ke.8, 12 = Ke.4, 16 = -0.006628571
Ke.1, 8 = 135โยทโ(11โยทโฮฑยฒ โ 132โยทโฮฑยฒ) + 5โยทโฮฝ + 150 = 135โยทโ(11โยทโ1 โ 132โยทโ1) + 5โยทโ0.2 + 150 = 0.1685714
Ke.2, 7 = -Ke.1, 8 = -0.1685714 , Ke.3, 6 = Ke.1, 8 = 0.1685714 , Ke.4, 5 = -Ke.1, 8 = -0.1685714 , Ke.9, 16 = Ke.1, 8 = 0.1685714
Ke.10, 15 = -Ke.1, 8 = -0.1685714 , Ke.11, 14 = Ke.1, 8 = 0.1685714 , Ke.12, 13 = -Ke.1, 8 = -0.1685714
Ke.1, 9 = -โ5435โยทโ(ฮฑa12 + 1ฮฑโยทโb12) + 7225โยทโA1 = -โ5435โยทโ(10.62 + 11โยทโ0.62) + 7225โยทโ0.36 = -0.5714286
Ke.5, 13 = Ke.1, 9 = -0.5714286
Ke.1, 10 = 135โยทโ(27โยทโฮฑa1 + 13ฮฑยฒโยทโb1) โ 625โยทโb1 = 135โยทโ(27โยทโ10.6 + 131โยทโ0.6) โ 625โยทโ0.6 = 1.504762
Ke.2, 9 = -Ke.1, 10 = -1.504762 , Ke.5, 14 = -Ke.1, 10 = -1.504762 , Ke.6, 13 = Ke.1, 10 = 1.504762
Ke.1, 11 = 135โยทโ(13โยทโฮฑยฒa1 + 27ฮฑโยทโb1) โ 625โยทโa1 = 135โยทโ(13โยทโ10.6 + 271โยทโ0.6) โ 625โยทโ0.6 = 1.504762
Ke.3, 9 = -Ke.1, 11 = -1.504762 , Ke.5, 15 = Ke.1, 11 = 1.504762 , Ke.7, 13 = -Ke.1, 11 = -1.504762
Ke.1, 12 = -โ1370โยทโ(ฮฑยฒ + 1ฮฑยฒ) + 150 = -โ1370โยทโ(1 + 11) + 150 = -0.3514286
Ke.2, 11 = -Ke.1, 12 = 0.3514286 , Ke.3, 10 = -Ke.1, 12 = 0.3514286 , Ke.4, 9 = Ke.1, 12 = -0.3514286 , Ke.5, 16 = -Ke.1, 12 = 0.3514286
Ke.6, 15 = Ke.1, 12 = -0.3514286 , Ke.7, 14 = Ke.1, 12 = -0.3514286 , Ke.8, 13 = -Ke.1, 12 = 0.3514286
Ke.1, 13 = 235โยทโ(27โยทโฮฑa12 โ 78ฮฑโยทโb12) โ 7225โยทโA1 = 235โยทโ(27โยทโ10.62 โ 781โยทโ0.62) โ 7225โยทโ0.36 = -16.095238
Ke.5, 9 = Ke.1, 13 = -16.095238
Ke.2, 16 = 135โยทโ(2โยทโa1ฮฑยฒ โ 13โยทโฮฑโยทโb13) + 2โยทโa175 = 135โยทโ(2โยทโ0.61 โ 13โยทโ1โยทโ0.63) + 2โยทโ0.675 = -0.024
Ke.4, 14 = -Ke.2, 16 = 0.024 , Ke.6, 12 = Ke.2, 16 = -0.024 , Ke.8, 10 = -Ke.2, 16 = 0.024
Ke.3, 3 = 435โยทโ(ฮฑ3 + 13ฮฑ) + 8โยทโฮฑ25 = 435โยทโ(13 + 131) + 8โยทโ125 = 1.92
Ke.7, 7 = Ke.3, 3 = 1.92 , Ke.11, 11 = Ke.3, 3 = 1.92 , Ke.15, 15 = Ke.3, 3 = 1.92
Ke.3, 4 = 235โยทโ(ฮฑยฒโยทโb1 + 11โยทโa13โยทโฮฑ) + 2โยทโb1โยทโโ(โ5โยทโฮฝ + 1โ)โ75 = 235โยทโ(1โยทโ0.6 + 11โยทโ0.63โยทโ1) + 2โยทโ0.6โยทโโ(โ5โยทโ0.2 + 1โ)โ75 = 0.192
Ke.7, 8 = -Ke.3, 4 = -0.192 , Ke.11, 12 = -Ke.3, 4 = -0.192 , Ke.15, 16 = Ke.3, 4 = 0.192
Ke.3, 7 = 235โยทโ(-2โยทโฮฑ3 + 9ฮฑ) โ 8โยทโฮฑ25 = 235โยทโ(-2โยทโ13 + 91) โ 8โยทโ125 = 0.08
Ke.11, 15 = Ke.3, 7 = 0.08
Ke.3, 8 = 135โยทโ(2โยทโฮฑยฒโยทโb1 โ 13โยทโa13โยทโฮฑ) + 2โยทโb175 = 135โยทโ(2โยทโ1โยทโ0.6 โ 13โยทโ0.63โยทโ1) + 2โยทโ0.675 = -0.024
Ke.4, 7 = -Ke.3, 8 = 0.024 , Ke.11, 16 = -Ke.3, 8 = 0.024 , Ke.12, 15 = Ke.3, 8 = -0.024
Ke.3, 11 = 335โยทโ(ฮฑ3 + 3ฮฑ) + 2โยทโฮฑ25 = 335โยทโ(13 + 31) + 2โยทโ125 = 0.4228571
Ke.7, 15 = Ke.3, 11 = 0.4228571
Ke.3, 12 = -โ170โยทโ(3โยทโฮฑยฒโยทโb1 + 13โยทโa13โยทโฮฑ) โ b1150 = -โ170โยทโ(3โยทโ1โยทโ0.6 + 13โยทโ0.63โยทโ1) โ 0.6150 = -0.06685714
Ke.4, 11 = -Ke.3, 12 = 0.06685714 , Ke.7, 16 = -Ke.3, 12 = 0.06685714 , Ke.8, 15 = Ke.3, 12 = -0.06685714
Ke.3, 15 = 135โยทโ(26ฮฑ โ 3โยทโฮฑ3) โ 2โยทโฮฑ25 = 135โยทโ(261 โ 3โยทโ13) โ 2โยทโ125 = 0.5771429
Ke.7, 11 = Ke.3, 15 = 0.5771429
Element stiffness matrix coefficients (above the main diagonal only)
Ke = D1,1โยทโKe = 24.305556โยทโKe = 796.296296135.185185135.18518516.736111-391.20370484.953704-13.6574074.097222-13.88888936.57407436.574074-8.541667-391.203704-13.65740784.9537044.097222 046.66666721.5972224.666667-84.95370414.027778-4.0972220.7083333-36.57407410.2777788.541667-1.625-13.6574071.9444444.097222-0.5833333 0046.6666674.666667-13.6574074.0972221.944444-0.5833333-36.5740748.54166710.277778-1.625-84.953704-4.09722214.0277780.7083333 0000.9777778-4.0972220.70833330.5833333-0.1611111-8.5416671.6251.625-0.2305556-4.0972220.58333330.7083333-0.1611111 0000796.296296-135.185185135.185185-16.736111-391.20370413.65740784.953704-4.097222-13.888889-36.57407436.5740748.541667 0000046.666667-21.5972224.66666713.6574071.944444-4.097222-0.583333336.57407410.277778-8.541667-1.625 00000046.666667-4.666667-84.9537044.09722214.027778-0.7083333-36.574074-8.54166710.2777781.625 00000000.97777784.0972220.5833333-0.7083333-0.16111118.5416671.625-1.625-0.2305556 00000000796.296296-135.185185-135.18518516.736111-391.203704-84.95370413.6574074.097222 00000000046.66666721.597222-4.66666784.95370414.027778-4.097222-0.7083333 000000000046.666667-4.66666713.6574074.0972221.9444440.5833333 000000000000.9777778-4.097222-0.7083333-0.5833333-0.1611111 000000000000796.296296135.185185-135.185185-16.736111 000000000000046.666667-21.597222-4.666667 0000000000000046.6666674.666667 0000000000000000.9777778
Element load vector
โFe = qโยทโA124โยทโ[6; a1; b1; A16; 6; -a1; b1; -A16; 6; -a1; -b1; A16; 6; a1; -b1; -A16] = 10โยทโ0.3624โยทโ[6; 0.6; 0.6; 0.366; 6; -0.6; 0.6; -0.366; 6; -0.6; -0.6; 0.366; 6; 0.6; -0.6; -0.366] = [0.9โ0.09โ0.09โ0.009โ0.9โ-0.09โ0.09โ-0.009โ0.9โ-0.09โ-0.09โ0.009โ0.9โ0.09โ-0.09โ-0.009] kN
Global stiffness matrix
K = 1020135.185185135.18518516.736111-391.203704-13.65740784.9537044.097222000000000000โฏ0 135.18518546.66666721.5972224.666667-13.6574071.9444444.097222-0.5833333000000000000โฏ0 135.18518521.59722246.6666674.666667-84.953704-4.09722214.0277780.7083333000000000000โฏ0 16.7361114.6666674.6666670.9777778-4.0972220.58333330.7083333-0.1611111000000000000โฏ0 -391.203704-13.657407-84.953704-4.0972221592.592593270.3703700-391.203704-13.65740784.9537044.09722200000000โฏ0 -13.6574071.944444-4.0972220.5833333270.3703793.33333300-13.6574071.9444444.097222-0.583333300000000โฏ0 84.9537044.09722214.0277780.70833330093.3333339.333333-84.953704-4.09722214.0277780.708333300000000โฏ0 4.097222-0.58333330.7083333-0.1611111009.3333331.955556-4.0972220.58333330.7083333-0.161111100000000โฏ0 0000-391.203704-13.657407-84.953704-4.0972221592.592593270.3703700-391.203704-13.65740784.9537044.0972220000โฏ0 0000-13.6574071.944444-4.0972220.5833333270.3703793.33333300-13.6574071.9444444.097222-0.58333330000โฏ0 000084.9537044.09722214.0277780.70833330093.3333339.333333-84.953704-4.09722214.0277780.70833330000โฏ0 00004.097222-0.58333330.7083333-0.1611111009.3333331.955556-4.0972220.58333330.7083333-0.16111110000โฏ0 00000000-391.203704-13.657407-84.953704-4.0972221592.592593270.3703700-391.203704-13.65740784.9537044.097222โฏ0 00000000-13.6574071.944444-4.0972220.5833333270.3703793.33333300-13.6574071.9444444.097222-0.5833333โฏ0 0000000084.9537044.09722214.0277780.70833330093.3333339.333333-84.953704-4.09722214.0277780.7083333โฏ0 000000004.097222-0.58333330.7083333-0.1611111009.3333331.955556-4.0972220.58333330.7083333-0.1611111โฏ0 000000000000-391.203704-13.657407-84.953704-4.0972221592.592593270.3703700โฏ0 000000000000-13.6574071.944444-4.0972220.5833333270.3703793.33333300โฏ0 00000000000084.9537044.09722214.0277780.70833330093.3333339.333333โฏ0 0000000000004.097222-0.58333330.7083333-0.1611111009.3333331.955556โฏ0 โฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฎโฑโฎ00000000000000000000โฏ0.9777778
Global load vector
โF = [0.9โ0.09โ0.09โ0.009โ1.8โ0.18โ0โ0โ1.8โ0.18โ0โ0โ1.8โ0.18โ0โ0โ1.8โ0.18โ0โ0โ...โ0.009] kN
Solution of the system of equations
โZ = slsolveโ(โK; โFโ)โ = [0โ0.5523352โ0.382696โ-0.4159901โ0.2028057โ0.3732616โ0.2648531โ-0.1936926โ0.2989084โ0.3091046โ0.04831304โ-0.02500182โ0.2612117โ0.3426307โ-0.1651507โ0.1275183โ0.1211517โ0.4681261โ-0.2671228โ0.2293407โ...โ-0.4159902] mm
Joint displacements
transpโ(โWzโ)โ = 00.3030.4880.5120.3830.16500.1390.340.4690.4720.3470.14600.1460.3470.4720.4690.340.139โฏ0 0.2030.4190.5620.5810.4850.3370.250.310.4380.5310.5330.4430.3110.2420.3110.4430.5330.5310.4380.31โฏ0.203 0.2990.4820.6050.620.5370.4170.350.3870.4850.5620.5640.4890.3870.3380.3870.4890.5640.5620.4850.387โฏ0.299 0.2610.4610.5930.6080.5130.3750.2990.3420.4550.5420.5440.460.3450.2890.3450.460.5440.5420.4550.342โฏ0.261 0.1210.3860.550.5650.4430.2530.1380.2170.3810.4930.4970.3890.2250.1340.2250.3890.4970.4930.3810.217โฏ0.121 00.340.5270.5420.4040.17300.1350.3360.4630.4670.3450.14500.1450.3450.4670.4630.3360.135โฏ0 0.1390.3980.5560.5660.440.2470.1290.2060.3670.4780.4810.3740.2110.1210.2110.3740.4810.4780.3670.206โฏ0.139 0.2990.4870.6080.6120.5110.3690.2870.3250.4330.5160.5190.4370.3240.270.3240.4370.5190.5160.4330.325โฏ0.299 0.3630.5260.6320.6350.5440.420.350.3760.4640.5360.5380.4670.3730.3290.3730.4670.5380.5360.4640.376โฏ0.363 0.2990.4870.6080.6120.5110.3690.2870.3250.4330.5160.5190.4370.3240.270.3240.4370.5190.5160.4330.325โฏ0.299 0.1390.3980.5560.5660.440.2470.1290.2060.3670.4780.4810.3740.2110.1210.2110.3740.4810.4780.3670.206โฏ0.139 00.340.5270.5420.4040.17300.1350.3360.4630.4670.3450.14500.1450.3450.4670.4630.3360.135โฏ0 0.1210.3860.550.5650.4430.2530.1380.2170.3810.4930.4970.3890.2250.1340.2250.3890.4970.4930.3810.217โฏ0.121 0.2610.4610.5930.6080.5130.3750.2990.3420.4550.5420.5440.460.3450.2890.3450.460.5440.5420.4550.342โฏ0.261 0.2990.4820.6050.620.5370.4170.350.3870.4850.5620.5640.4890.3870.3380.3870.4890.5640.5620.4850.387โฏ0.299 0.2030.4190.5620.5810.4850.3370.250.310.4380.5310.5330.4430.3110.2420.3110.4430.5330.5310.4380.31โฏ0.203 00.3030.4880.5120.3830.16500.1390.340.4690.4720.3470.14600.1460.3470.4720.4690.340.139โฏ0 mm
Bending moments
Zjโ(โjโ)โ = sliceโ(โโZ; k1โยทโโ(โj โ 1โ)โ + 1; k1โยทโjโ)โ
Zeโ(โeโ)โ = hpโ(โ[Zjโโ(โej.e, 1โ)โ; Zjโโ(โej.e, 2โ)โ; Zjโโ(โej.e, 3โ)โ; Zjโโ(โej.e, 4โ)โ]โ)โ
Average bending moments at joints, kNm/m
Mj = 1.5002750.30977470.21973790.15634340.15704780.99832250.15642330.15196360.19424650.1519630.15642380.99832240.15704870.15634320.21973730.30975011.5002588.502216.479335.777678โฏ1.500075 1.5613237.8059929.3381267.66833.333101-28.3619893.1463977.3238838.8950477.3238823.146398-28.3619893.3331027.66839.3381267.8059781.561340.32204035.3824847.126847โฏ1.561519 8.0886973.7662450.4861465-2.479522-4.4594030.15500784.7804022.842763-2.410353ร10-8-2.842763-4.780402-0.15500784.4594032.479522-0.4861463-3.766235-8.0886974.1090252.5724750.3669336โฏ8.088698
Bending moments for the plate
Bending moments - Mx
transpโ(โMxโ)โ = 1.5002758.5022111.06578310.4766316.7969380.8844119-31.2340010.30206295.6439378.7792398.8768585.919520.6860137-30.1316490.68601425.9195198.8768598.7792395.6439380.3020627โฏ1.500259 0.30977476.479339.4318288.8768314.922389-2.695046-10.242271-3.2600123.8085627.2499787.3609524.127543-2.764018-9.732553-2.7640184.1275437.3609527.2499783.808562-3.260011โฏ0.3097374 0.21973795.7776788.7161968.1349194.039416-2.257802-6.060553-2.8192132.932076.5165976.6403923.291629-2.294127-5.615407-2.2941273.2916296.6403926.5165962.93207-2.819213โฏ0.2197445 0.15634345.9931719.1259748.4927624.057756-3.257726-7.988419-3.8326612.9136446.7974156.9267753.304594-3.205142-7.381295-3.2051423.3045946.9267756.7974152.913644-3.832661โฏ0.1563406 0.15704787.2696610.372279.5566174.972566-5.268395-16.118845-5.8654883.7668247.7132127.8382354.157483-5.096247-15.085958-5.0962474.1574837.8382357.7132123.766824-5.865488โฏ0.1570494 0.99832259.03856711.08045310.1228085.763875-2.1353-38.650028-2.7455524.5217928.1908078.3150474.904107-2.088386-36.536389-2.0883864.9041078.3150478.1908074.521792-2.745552โฏ0.9983223 0.15642337.22106610.2820369.4364984.848272-5.357982-16.182847-5.9630023.624277.5601197.6988674.053292-5.138773-15.091479-5.1387734.0532927.6988677.5601193.62427-5.963002โฏ0.1564268 0.15196365.8732038.927978.2297863.773483-3.41909-8.032715-4.0090252.5942716.4697716.6264673.062933-3.269913-7.306549-3.2699133.0629336.6264676.4697712.594271-4.009025โฏ0.1519571 0.19424655.5141468.3802347.6929583.468366-2.491073-5.742609-3.073362.3104825.9792326.1426642.786926-2.376796-5.118358-2.3767962.7869266.1426645.9792322.310482-3.07336โฏ0.194254 0.1519635.8732048.927978.2297863.773483-3.41909-8.032715-4.0090252.5942716.4697716.6264673.062933-3.269913-7.306549-3.2699133.0629336.6264676.4697712.594271-4.009025โฏ0.1519564 0.15642387.22106610.2820369.4364984.848272-5.357982-16.182847-5.9630023.624277.5601197.6988674.053292-5.138773-15.091479-5.1387734.0532927.6988677.5601193.62427-5.963002โฏ0.1564266 0.99832249.03856711.08045310.1228085.763875-2.1353-38.650028-2.7455524.5217928.1908078.3150474.904107-2.088386-36.536389-2.0883864.9041078.3150478.1908074.521792-2.745552โฏ0.9983221 0.15704877.26965910.372279.5566174.972566-5.268395-16.118845-5.8654883.7668247.7132127.8382354.157483-5.096247-15.085958-5.0962474.1574837.8382357.7132123.766824-5.865488โฏ0.1570582 0.15634325.9931719.1259748.4927624.057756-3.257726-7.988419-3.8326612.9136446.7974156.9267753.304594-3.205142-7.381295-3.2051423.3045946.9267756.7974152.913644-3.832661โฏ0.156342 0.21973735.7776778.7161968.1349194.039416-2.257802-6.060553-2.8192132.932076.5165976.6403923.291629-2.294127-5.615407-2.2941273.2916296.6403926.5165962.93207-2.819213โฏ0.2197303 0.30975016.4793359.431838.8768314.922388-2.695046-10.242271-3.2600123.8085627.2499787.3609534.127543-2.764018-9.732553-2.7640194.1275437.3609527.2499783.808562-3.260012โฏ0.3094814 1.5002588.50222411.06578310.476636.7969390.8844112-31.2340010.3020635.6439378.7792398.8768585.919520.6860137-30.1316490.68601435.9195198.8768598.7792395.6439380.3020627โฏ1.500075 kNm/m
Bending moments My
transpโ(โMyโ)โ = 1.5613230.32204030.23543350.21672910.16938360.17865381.0129710.17816160.16785760.21141710.21080950.16727620.17645270.98644480.17645330.16727530.21081050.21141620.16785820.1781613โฏ1.561345 7.8059925.3824844.3092534.1357694.7426476.340088.2301146.282764.6121333.8928683.8680944.5307066.12227.9852796.12224.5307063.8680933.8928694.6121326.28276โฏ7.80597 9.3381267.1268475.8004925.5534226.4150978.115739.0465048.0146136.1889345.1451025.117036.0906867.821128.777037.821126.0906865.117035.1451016.1889348.014613โฏ9.338132 7.66835.4177814.2627213.9206234.2338295.4198786.3445995.3011323.9795653.478143.4810433.9686645.2228496.181445.2228493.9686643.4810433.478143.9795655.301132โฏ7.668298 3.333101-0.32799450.53974810.4057356-1.174273-3.412893-0.3301173-3.53113-1.3994770.099760860.1610125-1.210479-3.234594-0.2183684-3.234594-1.2104790.16101250.09976085-1.399477-3.53113โฏ3.333102 -28.361989-7.038398-1.936124-1.64658-4.99389-13.644591-36.317421-13.771031-5.215778-1.858062-1.763919-4.886758-13.033067-34.570998-13.033067-4.886758-1.763919-1.858062-5.215778-13.771031โฏ-28.361989 3.146397-0.50668750.34853730.1799829-1.45618-3.762354-0.7232665-3.922882-1.765717-0.2519029-0.1967723-1.596231-3.661551-0.6660923-3.661551-1.596231-0.1967723-0.2519029-1.765717-3.922882โฏ3.1464 7.3238835.0945673.9093.4891813.6755264.6859375.4988174.4840593.2542962.7944772.7886533.2145984.350075.2437544.350073.2145982.7886532.7944773.2542964.484059โฏ7.323878 8.8950476.7401375.3508914.9612985.6089926.9131347.5726116.6932525.1418814.1806234.1497385.0278956.4569637.2130836.4569635.0278954.1497384.1806235.1418816.693252โฏ8.895053 7.3238825.0945673.9089993.4891813.6755264.6859375.4988174.4840593.2542962.7944772.7886533.2145984.350075.2437544.350073.2145982.7886532.7944773.2542964.484059โฏ7.323877 3.146398-0.50668770.34853730.1799829-1.45618-3.762354-0.7232665-3.922882-1.765717-0.2519029-0.1967723-1.596231-3.661551-0.6660923-3.661551-1.596231-0.1967723-0.2519029-1.765717-3.922882โฏ3.1464 -28.361989-7.038398-1.936124-1.64658-4.99389-13.644591-36.317421-13.771031-5.215778-1.858062-1.763919-4.886758-13.033067-34.570998-13.033067-4.886758-1.763919-1.858062-5.215778-13.771031โฏ-28.36199 3.333102-0.32799460.53974820.4057355-1.174273-3.412893-0.3301173-3.53113-1.3994770.099760860.1610125-1.210479-3.234594-0.2183684-3.234594-1.2104790.16101250.09976085-1.399477-3.53113โฏ3.333109 7.66835.4177814.2627213.9206234.2338295.4198786.3445995.3011323.9795653.478143.4810433.9686645.2228496.181445.2228493.9686643.4810433.478143.9795655.301132โฏ7.6683 9.3381267.1268455.8004915.5534226.4150978.115739.0465048.0146136.1889345.1451025.117036.0906867.821128.777037.821126.0906865.117035.1451016.1889348.014613โฏ9.338126 7.8059785.3824794.3092544.135774.7426476.3400818.2301146.282764.6121333.8928683.8680944.5307066.12227.9852796.12224.5307063.8680933.8928694.6121326.28276โฏ7.805828 1.561340.32206640.23543390.21672770.16938450.17865291.0129710.17816180.16785750.21141720.21080940.16727630.17645270.98644480.17645340.1672750.21081080.2114160.16785810.1781614โฏ1.561519 kNm/m
Bending moments Mxy
transpโ(โMxyโ)โ = 8.0886974.1090251.416336-0.9453079-3.232331-4.7813530.021763664.8385873.3347481.141499-1.025238-3.194104-4.6732159.790543ร10-94.6732153.1941041.025238-1.141499-3.334748-4.838587โฏ-8.088697 3.7662452.5724750.9829256-0.515595-2.056503-3.2286130.058255113.3545242.2129550.7309473-0.6696087-2.124156-3.225638-2.241521ร10-83.2256382.1241560.6696087-0.7309473-2.212955-3.354524โฏ-3.766235 0.48614650.36693360.210340.08576266-0.1194038-0.31368240.07347950.46113970.26652590.05675809-0.07980332-0.2737226-0.4316644.761034ร10-90.4316640.27372260.07980332-0.05675809-0.2665259-0.4611397โฏ-0.4861468 -2.479522-1.796141-0.6013640.61671451.7469222.092690.07626484-1.948268-1.634146-0.58150580.49196681.5363311.879395-1.027996ร10-9-1.879395-1.536331-0.49196680.58150581.6341461.948268โฏ2.479522 -4.459403-3.240091-0.89783760.68261572.3709064.5250.07575709-4.382355-2.264704-0.68178660.60077812.1521384.2109132.395547ร10-10-4.210913-2.152138-0.60077810.68178662.2647044.382355โฏ4.459403 0.15500780.16110120.14918770.12455210.096356070.077948450.073859910.069867520.051235710.0210446-0.008412859-0.02659053-0.02344541-5.804372ร10-110.023445410.026590530.008412859-0.0210446-0.05123571-0.06986752โฏ-0.1550078 4.7804023.5659761.18935-0.4505739-2.206794-4.3993290.063467544.5354632.3793750.7252099-0.6264028-2.225368-4.2792951.444385ร10-114.2792952.2253680.6264028-0.7252099-2.379375-4.535463โฏ-4.780401 2.8427632.1343930.8682899-0.4377384-1.679619-2.0840940.037778392.1683581.7961380.6296579-0.5417111-1.677943-2.037693-3.548115ร10-122.0376931.6779430.5417111-0.6296579-1.796138-2.168358โฏ-2.842763 -2.410353ร10-83.062399ร10-8-4.397001ร10-93.066824ร10-10-8.12632ร10-113.636829ร10-11-1.325189ร10-114.764349ร10-12-7.14584ร10-131.790365ร10-13-3.8029ร10-132.910935ร10-1301.375317ร10-13-9.949168ร10-132.393ร10-12-3.712996ร10-122.380218ร10-123.027025ร10-12-1.394228ร10-11โฏ3.938371ร10-8 -2.842763-2.134393-0.86828990.43773841.6796192.084094-0.03777839-2.168358-1.796138-0.62965790.54171111.6779432.0376934.127238ร10-12-2.037693-1.677943-0.54171110.62965791.7961382.168358โฏ2.842763 -4.780402-3.565976-1.189350.45057392.2067944.399329-0.06346754-4.535463-2.379375-0.72520990.62640282.2253684.279295-1.911611ร10-11-4.279295-2.225368-0.62640280.72520992.3793754.535463โฏ4.780401 -0.1550078-0.1611013-0.1491877-0.1245521-0.09635607-0.07794845-0.07385991-0.06986752-0.05123571-0.02104460.0084128590.026590530.023445417.954439ร10-11-0.02344541-0.02659053-0.0084128590.02104460.051235710.06986752โฏ0.1550075 4.4594033.2400910.8978376-0.6826157-2.370906-4.525-0.075757094.3823552.2647040.6817866-0.6007781-2.152138-4.210913-3.335117ร10-104.2109132.1521380.6007781-0.6817866-2.264704-4.382355โฏ-4.459401 2.4795221.7961410.601364-0.6167145-1.746922-2.09269-0.076264841.9482681.6341460.5815058-0.4919668-1.536331-1.8793951.414276ร10-91.8793951.5363310.4919668-0.5815058-1.634146-1.948268โฏ-2.479516 -0.4861463-0.3669331-0.21034-0.085762740.11940380.3136824-0.07347951-0.4611397-0.2665259-0.056758090.079803320.27372260.431664-6.33653ร10-9-0.431664-0.2737226-0.079803320.056758090.26652590.4611397โฏ0.4861443 -3.766235-2.572475-0.98292630.51559492.0565033.228613-0.05825506-3.354524-2.212955-0.73094730.66960872.1241563.2256382.942834ร10-8-3.225638-2.124156-0.66960870.73094732.2129553.354524โฏ3.766119 -8.088697-4.109035-1.4163360.94530843.2323314.781353-0.0217637-4.838587-3.334748-1.1414991.0252383.1941044.673215-1.247931ร10-8-4.673215-3.194104-1.0252381.1414993.3347484.838587โฏ8.088698 kNm/m
@hydrostructai.com